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We propose a novel quasiparticle interpretation of the equation of state 
of deconfined QCD at finite temperature. Using appropriate thermal 
masses, we introduce a phenomenological parametrization of the onset 
of confinement in the vicinity of the phase transition. Lattice results 
^ of the energy density, the pressure and the interaction measure of pure 

SU(3) gauge theory are well reproduced. A relation between the ther- 
mal energy density of the Yang-Mills vacuum and the chromomagnetic 
q condensate (B 2 )t is found. We also present the two flavour QCD equa- 

CN ■ tion of state for realistic quark masses and apply the model to dilepton 

production in ultrarelativistic heavy-ion collisions. 

p ■ 1 Introduction 

^ . QCD undergoes a transition from a confined hadronic phase to a decon- 

fined partonic phase, the quark-gluon plasma (QGP), at a temperature of 
T c ~ 150 — 170 MeV [p]]. A central quantity of matter in thermal equilib- 
rium is the Helmholtz free energy from which the pressure p, energy density 
e and entropy density s - which are important ingredients for the description 
of ultra-relativistic heavy-ion collisions - can be derived. A first-principles un- 
derstanding of the equation of state (EOS) of hot QCD is therefore of great 
interest, also in order to reliably identify and calculate experimental signatures 
of that elusive state. Perturbative results on the EOS are available up to order 
0(g b s ) @. However, for temperatures of interest the strong coupling constant 
is presumably large: g s ~ 1 — 2. The perturbative expansion in powers of g s 
shows bad convergence already for much smaller values of the coupling. Fur- 
thermore, in the vicinity of a phase transition, perturbative methods are in 
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general not expected to be applicable. Non-perturbative methods such as lat- 
tice QCD calculations become mandatory. From these numerical simulations 
the EOS of a pure gluon plasma is known to high accuracy 0, and there are 
first estimates for the continuum EOS of systems including quarks, albeit still 
with unphysically large masses Q. 

Various interpretations of the lattice data have been attempted, most 
prominently as the EOS of a gas of quark and gluon quasiparticles. In a 
phenomenological framework, quarks and gluons are simply treated as non- 
interacting, massive quasiparticles ||. More recently, a quasiparticle descrip- 
tion of QCD thermodynamics has been derived in a more rigorous treatment 
using resummed hard thermal loop (HTL) perturbation theory ||. Employ- 
ing the full HTL spectral representions, the resulting EOS can be matched to 
lattice data down to temperatures T ~ 3 T c ; below that temperature, non- 
perturbative physics not amenable in an expansion in g s becomes important. 
Unfortunately, as evident from figure [3] (left panel), current experiments only 
probe that very temperature regime where the underlying physics, the confine- 
ment and chiral symmetry breaking mechanism, is not sufficiently understood. 
Phenomenological models incorporating as much physics as is known are there- 
fore necessary. Here, we propose a new quasiparticle model of the QGP that 
incorporates a parametrization of confinement close to T c , supplemented by 
thermal masses compatible with lattice results. Details can be found in 0. 

2 Quasiparticles and confinement 

The use of a quasiparticle model in QCD is based on the observation that in 
a strongly interacting system, the complex dynamics often rearranges itself in 
such a way that gross features of the physics can be described in terms of ap- 
propriate effective degrees of freedom. Consider a SU(3) gluon plasma. From 
asymptotic freedom, we expect that at very high temperatures the plasma 
consists of quasifree gluons. As long as the spectral function of the ther- 
mal excitations at lower temperatures resembles qualitatively this asymptotic 
form, a gluonic quasiparticle description is expected to be applicable. Their 
dispersion relation reads 

LU 2 k ^k 2 + ml(T), (1) 

where m*(T) acts as an effective mass generated dynamically by the interaction 
of the gluons with the heat bath background. 

However, the picture of a simple massive gas is presumably not appropri- 
ate close to T c because the driving force of the transition, the confinement 
process, is not taken into account. This physics has to be inserted by hand. 
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Figure 1: Left panel: Sketch of theory vs. lattice and the temperatures probed 
by current experiments. Right panel: Normalized e, s and p (solid lines) 
compared to continuum extrapolated SU(3) lattice data (symbols) ||. 

Below T c , the relevant degrees of freedom in SU(3) gauge theory are heavy, 
colour singlet glueballs. Approaching T c , deconfinement sets in and the gluons 
are liberated, followed by a sudden increase in entropy and energy density. 
Conversely, when approaching the phase transition from above, the number 
of thermally active degrees of freedom is reduced due to the onset of confine- 
ment. As T comes closer to T c , an increasing number of gluons gets trapped 
in glueballs which disappear from the thermal spectrum: since vtigb ~ 1.5 
GeV and T c ~ 270 MeV (for pure gauge theory), glueballs are simply too 
heavy to become thermally excited in the temperature range under considera- 
tion (up to about 5 T c ). So, while the confinement mechanism as such is still 
not understood, it is not necessary to know it in detail since we consider a 
statistical system. All confinement does on a large scale is to cut down the 
number of thermally active gluons as the temperature is lowered. This ef- 
fect can be included in the quasiparticle picture by modifying the distribution 
function of the gluons by a temperature-dependent confinement factor C(T): 

fBM^C(T)f B (LU k ). 

To become quantitative, we have to specify the thermal masses m*(T) en- 
tering (|l|) . Based on the observation that the Debye screening mass m d evalu- 
ated on the lattice shows approximate critical behaviour ||, in accordance with 
a weakly first order phase transition and in contrast to perturbative results, 
we parametrize 

m*(T) ~ G T (l - ^) . (2) 



Consider now the entropy of a gas of massive gluons with such a dropping 
m*(T). The result for s(T) will clearly overshoot the lattice entropy because 
light masses near T c lead to an increase in s(T). However, since the entropy is 
a measure for the number of active degrees of freedom, the difference may be 
accounted for by the aforementioned confinement process as it develops when 
the temperature is lowered towards T c . The explicit temperature dependence 
of the confinement factor C(T) can be obtained simply as the ratio of the 
lattice entropy and the entropy calculated with a dropping input gluon mass, 
and C(T) again shows near-critical behaviour: C(T) ~ (1 — T c /T) 7 . 

Thermodynamical quantities like the energy density can now be calculated as 

N 2 -l 7 

e(T) = \ dk k 2 [C(T)f B (u k )] u k + B{T). (3) 



The function B(T) is not an independent function, but uniquely determined 
by m* (T) , C (T) and their T-derivatives. It is necessary to maintain thermody- 
namical self-consistency and can be interpreted as the thermal energy density 
of the vacuum. The explicit expressions for B(T), the pressure and the en- 
tropy density can be found in [0. In figure [l] (right panel), we compare e, s 
and p with lattice data, as a function of T. We achieve a good and economic 
parametrization. The so-called interaction measure, related to the trace of 
the energy-momentum tensor, = e — 3p, is also nicely reproduced. Figure 
H (left panel) shows the function B{T) and the spacelike plaquette expecta- 
tion value A CT , as measured on the lattice, that can be related to the thermal 
chromomagnetic condensate (B 2 )t- 
We find the simple relation 

B(T) = l -K{T)T* = -^(B 2 ) T + V) T =o 

with the zero-temperature condensate (G 2 )t=o- This relation between B(T) 
and (B 2 )t may hint at a deeper connection between B(T) as a carrier of non- 
perturbative effects, and the magnetic condensate. After all, B(T) represents 
the thermal energy of the (non-trivial) Yang-Mills vacuum. Note that there is 
no such relation in other quasiparticle models, as evident from the dashed line 
in figure ^ (left panel). 

3 Dynamical quarks 

The extension of the mechanism presented so far to systems with dynamical 
quarks is not straightforward. Simulations of fermions on the lattice are still 




Figure 2: Left panel: The function B(T). Symbols show ~A CT |£§, the dashed 
line displays B(T) in other quasiparticle models ||. Right panel: Normalized 
e, s and p (solid lines) for Nf = 2 in the chiral limit. The grey band is a 
corresponding lattice estimation [10|. 



plagued by problems. However, when plotting the lattice pressure, normalized 
to the ideal gas value, for the pure gauge system and for systems with 2, 
2+1 and 3 quark flavours, it is found that the QCD EoS shows a remarkable 
flavour independence when plotted against T/T c . The flavour dependence is 
then well approximated by a term reminiscent of an ideal gas p(T, Nf) ~ 
(16 + 10.5 Nf)p(T/T c ) with a universal function p(T/T c ) ||. This hints at a 
confinement mechanism being only weakly flavour- dependent, and hence we 
assume that the function C(T) acts in a universal way on quarks and gluons. 
Figure |2] (right panel) shows pressure, energy and entropy density for two 
massless quark flavours. Reassuringly, the pressure of the confinement model 
lies well within a narrow lattice estimate for the continuum EOS in the chiral 
limit [[[(J. The picture for two light and a heavier strange quark looks similar 
although the approach to the ideal gas limit is slower. 

The neglect of the confined particles now violates detailed balance more 
severely since the lightest Goldstone modes have masses comparable to T c , 
hence they do contribute to the thermodynamics. Within a statistical hadroniza- 
tion model it may be possible to estimate their contribution. Furthermore, the 
calculation of susceptibilities will enable us to gain further insight in the nature 
of C(T). Work along these lines is in progress. 



4 Dilepton rates 



We apply the proposed model now to dilepton production in URHIC where it 
enters in two ways: the EOS serves as input for the construction of the fireball 



that is discussed in [II]]. Second, the dilepton emissivity of a static hot spot, 

dMee p v (u,k;T) 
d A xdujd?k ~ exp(cj/T) - 1' U 

is proportional to the vector spectral function py that is calculated from the 
photon self energy. By construction, the quasiparticle qq-loop is the only 
contribution. In addition, the confinement mechanism reduces the probability 
of thermally exciting a quark by C(T) and thus the total rate by C(T) 2 . The 



corresponding p v agrees with qualitative lattice results |L2 . 

To compare with experiment, the rate has to be convoluted with the 
fireball expansion. A discussion of the fireball dynamics and the final dilepton 



rates can be found in T. Renk's proceedings contribution fTT 
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